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fixed  narrow  range  of  encounter-pair  geometries  contribute  to  kgjj..  While  the  form 
of  the  plots  are  not  greatly  different,  the  latter  more  sophisticated  treatment 
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electrochemical-exchange  reactions  are  also  Included. 
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ABSTRACT 

The  predicted  dependence  of  the  blmolecular  rate  constant  for  outer- 
sphere  electron  exchange,  upon  the  longitudinal  relaxation  time, 

for  Debye  solvents  is  examined  numerically  on  the  basis  of  a  suitably 
combined  rate  formulation  in  order  to  examine  the  manner  and  extent  to 
which  the  rate-solvent  friction  dependence,  of  particular  experimental 
significance,  should  be  sensitive  to  the  degree  of  donor-acceptor 
electronic  coupling  and  related  factors.  The  treatment  accounts  for  the 
contributions  to  from  a  spatial  distribution  of  reactant  pairs  as  well 
as  for  the  effects  of  donor -acceptor  interactions  upon  the  unimolecular 
rate  constant,  k^^  (s*^),  for  each  encounter  geometry.  The  latter  include 
the  influence  of  electronic  interactions  as  prescribed  by  the  matrix 
coupling  element,  Hj2.  upon  the  effective  frequency  for  adiabatic  barrier 
crossing,  as  well  as  upon  the  electronic  transmission  coefficient, 

(i.e.  the  degree  of  reaction  nonadiabaticity) .  The  anticipated  dependence 
of  the  free-energy  barrier,  as  well  as  and  upon  the  donor-acceptor 

separation  is  accounted  for  in  the  k^^  calculations.  The  form  of  the 

resulting  log  k^^^  -  log  r'^  plots  exhibit,  as  expected,  a  marked  dependence 
on  the  value  of  for  reactant  contact,  Over  the  range 

appropriate  for  common  solvents  at  ambient  temperatures,  ca  1  x  10^^  to 
5  X  10^^  s'^,  the  log  k^j^  -  log  slopes,  x,  are  significantly  (or 

substantially)  below  unity  for  Hjj  <  0.5  kcal  mol"^;  x  decreases  markedly 

with  increasing  and/or  for  decreasing  .  The  approach  to  the  solvent 
inertial  limit  (for  large  r'^)  also  tends  to  decrease  x  further.  Numerical 
comparisons  are  made  with  corresponding  log  k^^  -  log  plots  obtained 
using  the  simple  "encounter  preequilibrium"  treatment  which  presumes  that 
only  a  fixed  narrow  range  of  encounter-pair  geometries  contribute  to  . 
While  the  form  of  the  plots  are  not  greatly  different,  the  latter  more 
sophisticated  treatment  generally  yields  smaller  slopes.  Some 


corresponding  numerical  calculations  for  electrochemical-exchange  reactions 
are  also  included. 


There  are  extensive  developments  occurring  at  the  present  time  in  our 
understanding  of  solvent  dynamical  factors  in  electron  transfer,  as  for 
other  condensed-phase  processes.^  The  original  treatments  for  electron 
transfer^  predict  that  the  barrier-crossing  frequency  for  adiabatic 
processes,  should  simply  be  proportional  to  the  inverse  of  the 

longitudinal  solvent  relaxation  time,  Subsequent  theoretical 

treatments  have  included  the  effects  of  barrier- top  curvature,^  non-Debye 
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solvent  relaxation,  and  reactant  vibrational  distortions,®  all  within 

the  dielectric  continuum  limit,  as  well  as  solvent  molecularity  effects  by 

employing  the  mean  spherical  approximation.® 

Most  experimental  kinetic  investigations  of  these  effects  involve 

examining  the  dependence  of  the  rate  parameters  for  bimolecular 

homogeneous -phase  or  electrochemical  electron-exchange  reactions  on  the 

solvent,  a  range  of  media  being  chosen  so  to  provide  a  suitably  large 

variation  in  the  solvent  dynamics . After  appropriate  correction  for 

the  solvent-dependent  barrier  height,  the  dependence  of  the  electron- 

exchange  rate  constant,  k^^^ ,  (or  the  inferred  barrier-crossing  frequency, 

upon  or  related  solvent  dynamical  parameters  can  yield  substantial 

insight  into  the  role  of  solvent  friction  upon  the  kinetics.  In  order 

for  such  solvent-friction  effects  to  be  manifested,  however,  it  is 

necessary  for  the  processes  to  remain  essentially  adiabatic  (or  nearly  so) 

within  the  range  of  nuclear  frequencies  accessed  by  the  solvent  relaxation 

dynamics.  We  recently  presented  numerical  calculations  of  the  dependence 

of  i/jj  upon  for  outer- sphere  electron-exchange  processes  having  varying 

degrees  of  donor-acceptor  orbital  overlap  as  described  by  the  electronic 

12 

coupling  matrix  element,  Hjj,  For  simplicity,  the  calculations  assumed 


2 


Debye  overdamped  relaxation,  the  effects  of  reaction  nonadiabaticity  being 

accounted  for  on  the  basis  of  a  generalized  Landau- Zener  treatment.  The 

results  indicate  that  substantial  values,  >  0.1  kcal  mol"^,  are 

required  to  maintain  reaction  adiabaticity  and  hence  preserve  the  full 

effect  of  solvent  friction,  as  Is  increased  within  the  range  ca  3  x 

10^°  to  3  X  10^^  s‘^  expected  for  common  polar  solvents  at  ambient 

12 

temperatures . 

The  calculations  presented  in  ref.  12  were  concerned  primarily  with 

the  solvent  frictional  dependence  of  the  "local"  unimolecular  rate 

constant,  ,  referring  to  a  particular  geometry  of  the  donor-acceptor 

pair.  Such  results  can  be  transposed  to  yield  corresponding  information 

for  experimental  bimolecular  (or  electrochemical)  rate  constants,  k^^^ ,  most 

13 

simply  by  utilizing  the  so-called  "encounter  preequilibrium"  treatment. 
This  approach  presumes  that  electron  transfer  takes  place  over  a 
sufficiently  narrow  range  of  precursor-complex  geometries  such  that  k^^^  can 
be  identified  approximately  with  k^^^  multiplied  by  a  "precursor  stability 
constant"  Kp.  However,  the  encounter  preequilibrium  treatment  is 
inherently  only  an  approximation  to  a  more  rigorous  description  of  k^^ 

involving  an  integration  of  k^^  over  a  distribution  of  reactant  separation 

14  15 

distances,  r.  As  alluded  to  recently  by  Beretan  and  Onuchic,  and  as 

outlined  in  detailed  fashion  herein,  the  relative  r-dependent  k^^ 

contributions  to  k^^^  are  inexorably  altered  somewhat  as  the  solvent 

friction  is  modified  (i.e.  is  varied)  for  a  given  reaction.  As  a 

consequence,  the  predicted  dependence  of  k^^^  upon  deduced  using  this 

spatial  integration  approach  may  differ  significantly  from  that  obtained 

from  the  simple  encounter  preequilibrium  formalism  that  utilizes  a  constant 

(i.e.  friction-independent)  Kp  term. 


3 


Given  the  central  importance  of  solvent -dependent  friction  effects 
upon  to  the  interpretation  of  the  experimental  kinetics ,  it  is 

crucial  to  elucidate  how  the  effects  of  electronic  coupling  combine  with 
solvent  friction  to  determine  k^^^  and,  especially,  its  dependence  upon  . 
The  objective  of  the  present  report  is  to  outline  in  detail  a  unified 
treatment  of  such  combined  electronic  coupling- solvent  friction  effects 
within  the  Debye  continuum  approximation  that  is  applicable  directly  to 
such  experimental  systems,  and  to  illustrate  its  application  by  means  of 
representative  numerical  calculations.  While  attention  is  focused 
primarily  on  homogeneous -phase  reactions,  some  pertinent  results  are 
included  for  electrochemical  exchange  processes.  We  present  in  a  companion 
paper  a  substantive  experimental  application  of  the  present  solvent 
friction-dependent  analysis,  so  to  provide  a  novel  means  of  extracting  Hj2 

Hr 

values  for  simple  self -exchange  reactions.  ° 

THEORETICAL  TREATMENT 
Bimolecular  Rate  Formulations 

As  noted  above,  the  overall  measured  rate  constant  for  bimolecular 

electron-exchange  reactions  in  homogeneous  solution,  k^^^ ,  strictly  is 

composed  of  an  integral  of  "local"  unimolecular  rate  constants,  kjt(i')- 
14 

Thus  we  can  write 

47rN  « 

*0 

where  N  is  Avogadro's  number,  r^  is  the  donor-acceptor  distance 
corresponding  to  the  (spherical)  reactant's  closest  approach,  and  gjjy^(r)  is 
the  radial  pair  distribution  function.  For  simplicity,  is  taken 


here  as  unity.  The  relation  corresponding  to  Eq.  (1)  that  describes 


electrochemical  rate  constants,  is 


-  S  k,t(r)  gp^(r)  dr 


(2) 


reflecting  the  linear  rather  than  spherical  spatial  coordinates  relevant  to 

V.  13c 

heterogeneous  process. 


The  local  rate  constant  generally  be  expressed  as 
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Kt  "  exp(-^AG  ) 


(3) 


where  -  (k^T)'  and  AG  is  the  free  energy  of  activation.  The  electronic 
transmission  coefficient  describes  the  extent  to  which  the  effective 
barrier-crossing  frequency  falls  below  the  nuclear  frequency  factor  as  a 
consequence  of  reaction  nonadiabaticity .  The  effect  of  solvent  friction  is 
manifested  entirely  in  the  term,  in  that  increases  as  the  solvent 
friction  (most  simply  r^)  decreases. 

All  three  components  of  k^^^  in  Eq.  (3):  and  AG  ;  are  expected 
to  depend  upon  the  reactant  separation  distance  r.  The  electronic  matrix 
coupling  element  will  decrease  markedly  with  increasing  r,  which  is 


usually  approximated  by  the  exponential  form: 
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-  (H°  exp[-Q(r  -  r.)] 


(A) 


where  Hjj  is  the  matrix  coupling  element  corresponding  to  r  -  r^ .  The 
coefficient  a  is  usually  calculated  (or  measured)  to  be  in  the  range  ca  1.0 
to  1.5 


The  dependence  of  upon  r  depends  on  the  magnitude  of  H°2  as  well 
as  that  of  a.  For  reactions  where  Hfj  is  sufficiently  large  so  that 
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reaction  adiabaticity  is  achieved,  -  1  for  r  -  r^^ .  As  r  is  increased 

beyond  r^,  however,  must  fall  eventually  to  the  point  where  electron 

transfer  becomes  nonadiabatlc  (l.e.  <  1),  whereby  the  Fermi  Golden  rule 

applies,  i.e.  a  (vide  infra). Under  such  nonadiabatic 

conditions,  however,  will  also  depend  upon  t/^;  for  «  1,  eventually 

a  so  that  the  net  preexponential  factor  becomes  independent 

12 

of  variations  in  the  solvent  friction. 

The  decrease  in  with  increasing  r  will  also  obviously  influence 

Ar*  • 

AG  since 

AG*  -  AG*  -  (5) 

where  AG^  is  the  barrier  for  "*  0.  i-e.  when  the  barrier  top  resembles  a 

cusp.  In  addition,  AG  is  expected  to  increase  sharply  with  increasing  r 

even  in  the  absence  of  such  "resonance  splitting"  since  according  to  the 

20 

dielectric  continuum  treatment  of  outer-shell  barriers: 


AG*-  (e2/4)(e;J  -  e'^)(a'^  -  r‘^) 


(6) 


where  e  is  the  electronic  charge,  and  are  the  optical  and  static 

(zero- frequency)  dielectric  constants,  and  a  is  the  reactant  radius.  The 

functional  dependence  of  AG*  upon  r'^,  predicted  by  Eq .  (6),  is  in  accord 

with  experimental  energies  of  tjptical  electron  transfer  within  mixed- 

21 

v.  lence  complexes  having  varying  internuclear  separations  r. 

The  relation  coiresponding  to  Eq.  (6)  for  electrochemical  reactions 


is 
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AG*  -  (eV8)(e;J  -  <;')((a-'  -  (2r,)-i] 


(6a) 
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where  is  the  distance  between  the  reactant  and  the  surface  plane.  Since 

outer-sphere  processes  are  often  anticipated  to  occur  at  the  outer 

Helmholtz  plane  where  the  reactant  is  separated  from  the  metal  surface  by  a 

solvent  monolayer,  the  influence  of  the  "imaging"  term  r^  upon  could  be 

unimportant , ‘  yielding  only  mild  increases  in  AG*  with  increasing  r^  . 

In  the  presence  of  solvent  friction,  is  anticipated  to  depend  on 
2-4 

the  barrier  height  as  well  as  on  the  shape  of  the  barrier  top,  and  hence 

2-4  3  12  22 

upon  Hj2-  Although  both  these  dependencies  are  relatively  mild,  ’  ’ 

since  tends  to  increase  both  with  increasing  AG^  and  decreasing  one 

can  predict  that  will  Increase  monotonically  with  increasing  r.  In  the 

transition-state  theory  (TST)  limit  (i.e.  in  the  absence  of  solvent 

friction),  however,  both  these  lependencies  will  disappear.  In  any  case, 

the  influence  of  these  factors  upon  the  h^^^-r  profile  are  overshadowed  by 

the  corresponding  increases  in  AG  and  decreases  in  with  increasing  r, 

so  that  k^j^(r)  will  inexorably  decrease  sharply  with  increasing  reactant 

separation . 

Provided  that  declines  sufficiently  rapidly,  to  a  first 

approximation  the  integral  in  Eq.  (1)  can  be  reexpressed  simply  in  the 

13 

"encounter  preequilibrium"  formulation  as 

\  Kt  -  \  *.i  ‘'n  exp(*^4G*)  (7) 

The  preequilibrium  constant  reflects  the  spatial  range  over  which 

significant  contributions  to  occur,  with  k^^  and  its  component 

parameters  and  AG*  being  approximated  by  the  local  values  at 

r  -  r^.  However,  as  noted  above  it  is  preferable  for  the  present  purposes 
to  perform  a  spatial  k^^^  integration  to  find  k^^^  since  Kp  is  expected  to 
change  somewhat  as  the  solvent  friction  is  altered  (vide  infra) . 


In  the  sections  that  follow,  we  outline  a  simple  combined  treatment  of 
the  above  effects  upon  followed  by  illustrative  numerical 
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predictions  of  the  resulting  dependence  of  upon  solvent  friction. 

The  Unimolecular  Rate  Constant 

We  require  here  a  treatment  of  which  accounts  for  the 

alterations  in  needed  for  the  complete  spatial  integration  in 

Eq.  (1),  together  with  the  ability  to  consider  wide  variations  in  the 

extent  of  solvent  friction.  Considered  first  is  a  suitable  means  of 

describing  barrier  shape  effects  upon  followed  by  a  formalism  which  can 

account  additionally  for  nonadiabaticity . 

In  the  presence  of  solvent  friction,  as  increases  and  the  barrier 

top  becomes  increasingly  rounded  the  frequency  of  successful  barrier 

crossings  is  progressively  diminished.  This  effect,  first  considered  by 
23 

Kramers,  has  also  been  described  recently  in  the  context  of  electron- 

3  ^  12 

transfer  processes.  ’  ’  In  our  previous  examination  of  -  solvent 

12 

friction  dependencies,  we  employed  the  Smoluchowski  treatment  described 

by  Calef  and  Wolynes .  This  approach  describes  the  passage  across  the 

potential-energy  surface  in  terms  of  simple  diffusion.  More  generally,  the 

diffusion  coefficient,  D,  can  be  considered  to  depend  on  the  position,  q, 

24  25 

along  the  reaction  coordinate  potential-energy  surface,  V(q).  ’  Whereas 

D  will  be  independent  of  q  for  a  pair  of  harmonic  surfaces  forming  a  cusp 

barrier,  when  >  0  the  diffusion  coefficient  within  the  intersection 

region,  will  differ  from  that  within  the  parabolic  wells,  D^.  Indeed, 

for  large  values  (i.e.  for  large  resonance  splittings),  the 

24 

component  dominates  the  overall  barrier-crossing  rate. 


An  expression  for  the  adiabatic  unlnolecular  rate  constant,  ,  due 

24 


to  Okuyana  and  Oxtoby  is 


-  2  /  dq[D(q)  P,_(q)]-M;  dq'  P.,<q')] 


(8) 


where  “  2  exp[-^  V(q)]  (8a) 

00 

and  Z  -  J*  exp(-^  V(q)]  dq  «  (8b) 


(The  approximation  in  Eq.  (8b)  arises  since  the  greatest  contribution  to 
the  integral  will  be  close  to  the  Initial  state  where  V(q)  -  0.5  m  q^) 
Tne  value  of  q  is  taken  to  be  zero  for  the  reactant  initial  state;  q 
refers  to  the  transition  state;  n  is  the  "effective  mass"  of  the  reacting 
particle,  and  is  the  well  frequency.  The  above  integrals  may  also  be 

simplified  for  large  barriers  (^AG*  >  5)  to  yield: 

« 

(k.'t)'"  -  2  /  dq[D(q)  P,,(q)]-'  (9) 


Equation  (9)  is  adequate  for  the  numerical  calculations  presented  below, 

*  23 

where  ^AG  ~  10.  This  relation  reduces  to  the  simple  Kramers  relation 

2  ^ 

(e.g.  ref.  3)  if  D  is  independent  of  q;  it  also  can  be  related  directly 
to  the  well-known  Grote -Hynes  result. 

In  order  to  apply  this  treatment  to  solvent  friction  in  electron 
transfer,  we  require  appropriate  expressions  for  both  V(q)  and  D(q).  The 
former  can  be  expressed 

V(q)  -  0.5[Vj(q)  +  V2(q)l  -  0.5{[V,(q)  -  V2(q)]2  -h  (10) 


where  Vj  and  Vj  refer  to  the  reactant  and  product  potential -energy  curves, 
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respectively.  These  quantities  can  be  expressed  in  terms  of  energetic  (E) 
as  well  as  position  (q)  coordinates.  Thus  Vj^(q)  -  0.5  m  q^,  whereas 
Vj^(E)  -  0.5  nij  E*,  where  -  (2Ej  Convenient  expressions  for  Vj 

and  in  the  latter  (energy)  coordinates  are 

V^(E)  -  eV^E,;  V^(E)  -  (E  -  2E^)V^E,  dD 

where  E^  is  the  reorganization  energy  (-  4AG  for  a  symmetric  barrier  as 
considered  here). 

When  ^  barrier  top  can  often  be  approximated  by  an  inverted 

parabola,  allowing  the  potential -energy  surface  close  to  the  transition 
region  to  be  expressed  as 

Vj^(E)  -  AG*  -  0.5  mj  ^1^2) 

where  the  "characteristic  barrier-top"  frequency  is  related  to  the  well 

22  26 

frequency  for  symmetric  barriers  by  ’ 

-  [(2AG*/Hi2)  -  1]^*  (13) 

Given  this  form  of  V(E),  a  general  expression  for  the  position- 
dependent  diffusion  coefficient  can  be  written  as  (see  Appendix  A): 

Dp(P  -  w,b)  -  (^  fflj  (14) 

where  is  the  effective  friction.  This  quantity  can  be  expressed  for 

Debye  solvents  as  (Appendix  A) : 

»?:“  -  2«f/q{l  -  (1  -  (Aw^V))**)  (for  V  >  2w„)  (15a) 


or 


»?:"  -  2w*/r, 


(for  q  <  2u^) 


(15b) 


-  2u>1/t,[[1  +  (4a)2/,,2)]^  -  1])  (16) 

where  t}  -  <*>*  r^.  Equations  (14) -(16)  enable  and  Dj^,  and  hence  to 

be  determined  in  terms  of  r^,  u^,  Hjj.  AG  which  together  describe  the 

experimental  reactant-solvent  system.  The  calculation  of  by  this  means 

employs  for  potential  energies  E  <  E,  and  for  E  >  E^ ,  where  E^ 

defines  the  energy  where  the  well  potential  surface  gives  way  to  the 

inverted  parabola  characteristic  of  the  barrier  top.  This  approach 

therefore  provides  a  relatively  general  as  well  as  convenient  means  of 

accounting  for  such  "barrier- top  roundedness",  caused  by  moderate  or  large 

electronic  coupling,  within  the  adiabatic  limit. 

As  noted  above,  the  effects  of  electronically  nonadiabatic  barrier 

crossing  also  need  to  be  accounted  for,  not  only  in  the  k^^^-r  integration 

to  find  k^jj  in  a  given  solvent  (Eq.  (1)]  but  also  so  to  provide  a  more 

complete  description  of  the  k^^^  -  dependence.  Whereas  kj^  is  predicted 

to  be  proportional  to  r*^  in  the  absence  of  solvent  inertial  effects  (i.e. 

for  overdamped  solvent  motion),^  k^^  will  fall  below  kj^^  for  small 

and/or  large  and  eventually  will  become  independent  of  .  One 

means  of  bridging  these  two  regimes  is  to  regard  the  overall  unimolecular 

reaction  to  be  composed  of  separable  rate  constants  for  approaching  and 

2fi  ^  3  l^fl 

then  crossing  the  barrier  top,  so  that  ’  ’  ’ 

(k„)-‘  -  (17) 

where  k^^  accounts  for  the  possible  influence  of  reaction  nonadiabaticity 
(i.e.  transitions  to  the  upper  diabatic  surface,  rather  than  adiabatic 
barrier  crossings  involving  electron  transfer).  The  effect  of  solvent 


friction  is  contained  wholly  within 

A  significant  limitation,  however,  of  Eq.  (17)  for  the  present 

purposes  is  that  it  is  strictly  valid  only  for  sufficiently  small  "barrier- 

top  roundedness"  (small  so  that  k*^^  is  dominated  by  motion  within  the 

potential-energy  wells  rather  than  in  the  intersection  region  (vide  infra). 

Equation  (17)  therefore  cannot  account  for  the  effects  of  barrier-top  shape 

upon  k“j^.  [Note  that  a  relation  of  apparently  similar  form  to  Eq.  (17)  has 

been  used  to  formulate  such  barrier-shape  effects  upon  k^^  within  the 

adiabatic  limit, ^  where  the  first  term  refers  to  motion  within  the  wells 

and  the  second  describes  the  barrier-crossing  step  (see  also  ref.  22).]  An 

expression  for  the  "barrier-crossing"  rate  constant,  k^^,  valid  at  least  in 

27 

the  normal  (as  opposed  to  inverted)  free-energy  regime  as  considered  here 
is 

“  ^ei^‘^o/2^)  exp(-/3AG*)  (18) 

where  denotes  the  rate  constant  in  the  TST  limit,  and  (not  to  be 

confused  with  the  overall  electronic  transmission  coefficient  employed 

2  8 

above)  can  be  expressed  as 

«0 

^.1  -  /  exp(-^E^)  d(;9E^)  (19) 

o 

where 

-  [1  -  exp(-27r7)  )/(l  -  0.5  exp(-2»r7)]  (19a) 

7  -  |H,JV2K(E,E,)‘‘  (19b) 

and  E^(-  E  -  AG*)  is  the  system  energy  relative  to  that  at  the  barrier 


12 


As  shown  in  Appendix  B,  Jc^j^  can  be  expressed  to  a  good  approximation 
more  simply  as: 


(20) 

V  19 

where 

7i  -  |H,JV2ft  w^(E,  kjT)** 

(20a) 

13a  14a 

Another,  more  familiar,  form  for  Jc^j^  is  ‘  (See  Appendix  B)  ; 


'‘•1  =  1  +  W, 


2(1  -  exp(-2^3/2^^jj 

2  -  exp(-2ff^^^7j) 


where  [  -  1  -  exp( -2jr^^^7j,)  ]  is  the  averaged  Landau- Zener  probability 
for  a  "single  crossing"  of  the  barrier  top.  Equation  (21)  is  obtained  as 
an  approximation  to  Eq.  (19)  by  replacing  the  7  term  in  [Eq.  (19a)]  by 
its  average  value  ?r^^^7j. 

From  Eqs .  (17)  and  (18)  we  can  express  the  distance -dependent 
unimolecular  rate  constant 


s 


1  +  [X.,  k^^Vk,?] 


(22) 


A  closely  related  simplified  expression  for  resulting  from  Eqs. 

(17)  and  (20),  utilized  previously, 


2it  |HjJ^(16ft2  m  AG*  exp(-)9AG*) 

1  +  (m  |HiJ2  r^/f,  AG*) 


Kt  - 


This  relation  is  derived  from  Eq.  (17)  by  assuming  that  k"^^  refers  to 

2-4 

overdamped  solvent  motion  over  a  cusp  barrier  (i.e.  -♦  0),  whereupon 


(/9AG*/4jr)^  exp(-/3AG*) 


(24) 
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and  «  kjT  so  that  the  expression  for  reduces  to 

k"J  -  2»r  exp(-^AG*)  (25) 

■  7^  k^®^  (25a) 

Although  Eq.  (23)  as  well  as  Eq.  (22)  provides  a  useful  approximate 

form  under  conditions  where  is  small,  it  does  not  account  either  for 

the  influence  of  the  barrier- top  curvature  upon  or  for  the  onset  of  the 

inertial  limit,  where  -*  w^/2ff  (vide  infra).  These  shortcomings  may  be 

avoided  by  using  in  place  of  Eq.  (17)  the  following  modified  interpolation 

30 

expression,  obtained  by  Straub  and  Berne  [Eq.  3.15  of  ref.  30a],  which 
can  be  written  as 


1  -  W 

+  (— • )  (k^”)-i  (26) 

By  assuming  the  relation  between  and  given  in  Eq.  (21),  and 

rearranging,  we  obtain 


k.t  - 


(27) 


Equation  (27)  differs  from  the  simplified  Eq.  (22)  by  the  presence  of 

an  additional  term  in  the  denominator.  The  two  expressions  will 

therefore  only  become  equivalent  either  when  Jc^j^  «  1  and/or  when 

k’®^/k“^  »  1;  i.e.  for  small  Hjj  and/or  high  friction.  Equation  (27)  is 

clearly  preferred  for  the  present  purposes  since  it  should  be  applicable 

lie 

also  to  systems  having  large  values,  a  situation  often  encountered  ° 
and  Indeed  required  in  order  for  solvent  dynamical  effects  to  influence  k^^^ 
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in  low  friction  media  (vide  supra).  A  related  point  is  that  Eq.  (27), 
unlike  Eq.  (22),  yields  the  correct  prefactor,  w^/2ir ,  in  the  adiabatic 
"zero-friction"  (i.e.  TST)  limit.  (Equation  (22)  yields  a  prefactor  that 
is  twofold  smaller  in  this  limit.] 

Equation  (27)  therefore  describes  in  a  straightforward  as  well  as 
general  fashion  the  separate  effects  of  reaction  nonadiabaticity  and 
solvent  friction  via  the  and  terms,  respectively.  The  influence  of 
solvent  inertia  appears  in  both  k^^^  and  k“^.  If  the  reaction  involves 
sufficient  electronic  coupling  to  be  entirely  adiabatic,  then  k^^^  -  k“^  in 
the  presence  of  solvent  friction,  or  k^^^  -  k^^^  in  its  absence.  Even 
though  Eq.  (27)  is  utilized  primarily  in  the  following  numerical 
calculations,  we  also  include  some  results  obtained  with  the  conventional 
simplified  formula  Eq .  (23)  for  comparative  purposes. 

The  form  of  the  preceding  parameters  can  readily  be  related  to  the 
phenomenological  quantities  /c^j^  and  in  Eq.  (3).  The  latter  parameter  is 
obtained  directly  from  [Eq.  (9)]  by  removing  the  Boltzmann  term  since 
[cf  Eq.  (3)] : 

Kx.  -  exp(-;9AG*)  (28) 

The  transmission  coefficient  can  be  identified  with  the  more 
familiar  quantity  in  the  TST  limit  where  solvent  friction  disappears 
such  that  2ni/^  “  “o’  presence  of  friction,  from  Eqs.  (3) 
and  (28)  can  be  expressed  simply  as  ^  general  relation 
between  and  can  be  obtained  simply  from  Eqs.  (21)  and  (25): 


(29) 


where  the  "adiabatic  transmission  coefficient"  «,  (  “  describes 

the  influence  of  solvent  friction  upon  in  the  absence  of  reaction 

12 

nonadiabaticity .  In  some  earlier  calculations,  we  accounted  for  the 

combined  effects  of  friction  and  nonadiabaticity  upon  k^^  by  using  a 

modified  form  of  Eq.  (21)  expressed  In  terms  of  rather  than  The 

numerical  similarities  of  these  treatments  are  noted  briefly  in  Appendix  B. 


NUMERICAL  RESULTS  AND  DISCUSSION 
Homogeneous  Bimolecular  Reactions 

Having  outlined  a  relatively  versatile  means  of  predicting  the 
functional  dependence,  at  least  in  the  continuum  Debye  limit,  it  is 
of  primary  interest  to  examine  the  extent  to  which  this  solvent- friction 
dependence  is  sensitive  to  the  fundamental  kinetic  parameters  involved,  and 
especially  to  the  extent  of  donor-acceptor  orbital  overlap,  H®2- 

Figure  1  displays  an  illustrative  set  of  log  k^^^-log  plots  for 
homogeneous  bimolecular  processes  obtained  using  the  conventional 
simplified  expression  for  k^^^  [Eq.  (23)]  along  with  Eqs.  (1),  (4)-(6),  for 
the  sequence  of  H^^  values  0.02,  0.05,  0.1,  0.2,  0.5,  and  1.0  kcal  mol'^. 
(The  kjjj  values  at  a  given  value  increase  with  increasing  Hj2-)  The 
other  parameters  required  in  this  analysis  were  held  constant  as  follows; 
AG*  -  5.0  kcal  mol"^,  o  -  1.25  A'^,  r^  -  2a  -  8  A  and  T  -  298  K.  These 
values  are  probably  typical  for  many  experimental  systems . 

Choosing  different  values  of  a  and  r^^,  at  least  within  the  ranges  normally 
expected  experimentally  (ca  0.8  to  1.5  A'^  and  3  to  6  A,  respectively)  have 
relatively  little  influence  on  the  form  of  these  plots.  Admittedly,  the 
larger  values  chosen  violate  the  assumption  Hj2  «  k^T  that  is  embodied 
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in  Eq.  (23);  It  Is  nevertheless  still  instructive  to  examine  the  influence 

of  altering  more  widely  upon  the  predicted  dependence. 

Figure  1  shows  that  while  approximately  a  at  relatively  small 

values,  as  expected,  increasing  t'^  leads  eventually  to  marked  decreases  in 

the  log  k^x*^°8  slope  and  eventually  to  a  solvent  friction- independent 

rate  in  accordance  with  the  Fermi  Golden  Rule,  whereupon  k^^^  q  • 

Necessarily,  then,  the  range  of  values  where  these  deviations  occur 

shift  markedly  upwards  as  H°2  is  increased;  i.e.  solvent  friction  effects 

persist  to  larger  values.  Even  for  and/or  values  that  are 

sufficiently  large  so  that  essentially  k^^^  a  ,  k^^^  is  seen  to  increase 

significantly  with  increasing  H°2 .  This  is  due  to  the  corresponding 

decreases  in  the  effective  barrier  height  [Eq.  (5)]. 

Figures  2A  6i  B  display  corresponding  sets  of  log  k^^^-log  plots 

obtained  by  utilizing  the  more  complete  k^„  expression  [Eq.  (27)], 

including  the  effects  of  barrier- top  roundedness  upon  These 

calculations  utilized  the  same  input  parameter  values  as  in  Fig.  1,  but 

with  inertial  limiting  (i.e.  TST)  frequencies,  ,  of  1  x  10^^  and 

5  X  10^^  s"^  in  Figs.  2A  and  B,  respectively.  (These  values  represent 

22 

typical  rough  lower  and  upper  limits  of  for  simple  polar  solvents.  ) 

As  before,  k^^^  was  calculated  from  k^t,(r)  using  Eqs .  (1),  (4) -(6),  with 

kj^(r)  obtained  using  Eqs.  (8) -(16),  from  Eq.  (20)  and  k^®^  from 

Eq.  (18).  The  evaluation  of  k^^^  In  this  manner  therefore  involves  a  double 

Integration,  both  spatially  [Eq.  (1)]  and  along  the  reaction  coordinate 

[Eq.  (8)].  The  double  integrals  were  evaluated  using  a  standard  two- 

31 

dimensional  quadrature  technique,  executed  on  the  Purdue  University  Cyber 
205  supercomputer  using  double  precision  (128  bit).  [The  single  integral 
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calculations,  for  example  in  Fig.  1,  were  undertaken  using  a  Zenith  387 

laboratory  microcomputer  with  6A  bit  precision  at  16  MHz.] 

The  overall  shapes  of  the  log  -  log  plots  in  Figs.  2A  and  B  do 

not  differ  greatly  from  those  in  Fig.  1,  indicating  that  the  inclusion  of 

the  effects  of  barrier- top  curvature  and  solvent  inertia  do  not  have  a 

dramatic  effect  on  the  rate-solvent  friction  dependence  under  most 

conditions.  Nevertheless,  the  influence  of  the  latte’"  factor  is  clearly 

seen  in  the  traces  in  Figs.  2A  and  B  for  the  hig>..' ?t  values  (i.e.  for 

the  strongest  electronic  coupling)  in  that  the  plots  deviate  from  linearity 

for  large  values  to  a  greater  extent  as  is  decreased  (compare  Fig. 

2A  versus  2B)  .  This  behavior  reflects  the  approach  to  the  inertial 

limiting  rate  as  the  solvent  friction  is  reduced  to  small  values  (i.e.  when 
32 

>  Wjj)  .  However,  in  some  experimental  systems  this  "inertial  limit" 

la  23 

may  not  be  observable,  ’  in  part  due  to  the  well-known  onset  of 

33 

diffusion  control.  The  primary  influence  of  barrier- top  roundedness  is 

to  diminish  the  displacements  between  the  log  k^^^  -  log  curves  in  the 
"adiabatic"  region,  corresponding  to  large  (i.e.  high  friction),  where 
roughly  k^^^  q  (compare  Figs.  2A,B  with  Fig.  1).  This  effect  arises 

because  the  inclusion  of  barrier-top  roundedness  (i.e.  deviations  from  a 
cusp)  tends  to  depress  and  therefore  k^^^  and  k^^^ ,  increasingly  as  H®^ 

and  hence  the  extent  of  this  curvature  increases  (vide  supra) . 

One  significant  feature  of  both  Figs.  1  and  2  is  that  the  slopes  of 
the  log  k^jj-log  plots  are  significantly  below  unity  even  within  the 
relatively  high  friction,  "adiabatic",  region.  Even  for  strong  electronic 
coupling,  say  when  Hjj  “  1  kcal  mol'^  (i.e.  the  uppermost  trace  in  these 
figures)  ,  the  log  k^^^  -  log  r'^  slopes  Increase  slightly  with  decreasing 
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For  example,  for  the  trace  where  H°2  -  1  kcal  mol'^  in  Fig.  2B,  the 
slope  increases  from  0.93  for  -  10^^  s'^  to  0.98  for  -  10®  s*^. 

Although  the  degree  of  nonlinearity  can  be  relatively  minor  (as  in  this 
example),  it  reflects  in  part  the  changes  in  the  form  of  the 
distribution  [i.e.  the  integrand  in  Eq.  (1)]  as  the  extent  of  solvent 
friction  is  altered  (cf  ref.  15). 

Given  this  effect,  it  is  of  interest  to  compare  more  closely  such 
log  kjjj  -  log  plots  obtained  by  performing  a  spatial  k^j_  integration 

[i.e.  using  Eq.  (1)]  with  corresponding  curves  resulting  from  the 
simplified  "encounter  preequilibrium"  treatment  [Eq.  (7)].  Figure  3 
contains  a  pair  of  comparisons  of  this  type  for  representative  Hj2  values 
of  0.05  and  0.5  kcal  mol'^.  In  both  cases,  the  solid  traces  were  extracted 
from  Fig.  2B,  i.e.  involving  a  spatial  k^^^  integration  and  including  the 
effects  of  both  barrier-top  roundedness  and  solvent  inertia.  The  dotted 
traces  were  taken  from  Fig.  1,  i.e.  with  only  spatial  k^^^  integration  taken 
into  account.  (The  numerical  parameters  used  in  Fig.  3  are  identical  to 
those  in  Figs.  1  and  2B.)  The  corresponding  dashed  curves  in  Fig.  3,  on 
the  other  hand,  include  the  effects  of  barrier- top  roundedness  and  solvent 
inertia  as  for  the  solid  traces,  but  k^^^  was  derived  from  k^^  by  employing 
Eq.  (7)  rather  than  using  the  spatial  integration  procedure.  The  k^^ 
values  used  to  calculate  k^^^  from  Eq.  (7)  refer  to  the  values  with  the 
reactants  essentially  in  contact,  i.e.  with  r  -  r^,  with  taken  as 
0.25  This  latter  value,  although  somewhat  arbitrary,  corresponds  to 

that  expected  for  ^  -  8  A,  with  an  effective  "reaction  zone 

thickness"  of  about  0.5  Since  is  presumed  to  be  independent  of 

both  and  Hjj  in  this  simple  analysis,  the  form  of  the  resulting 


log  -  log  curves  (dashed  traces  in  Fig.  3)  will  reflect  only  the 

friction- dependent  behavior  of  the  unimolecular  rate  constant  for  the 

reactants  in  contact  (cf  calculations  in  ref.  12). 

Although  in  some  respects  similar,  this  trio  of  plots  for  both 

H°2  -  0.05  and  0.5  kcal  mol"^  (lower  and  upper  sets,  respectively,  in 

Fig.  3)  show  subtle  yet  significant  differences.  In  particular,  the  curves 

calculated  using  the  spatial  k^^  integration  (solid,  dotted  traces)  exhibit 

significantly  greater  nonlinearities  as  well  as  smaller  slopes  than  that 

obtained  from  the  simple  preequilibrium  treatment.  Given  that  the 

values  for  simple  polar  solvents  at  ambient  temperatures  differ  by  only 

about  50  fold  (or  less),  within  the  range  ca  1  x  10^^  to  5  x  10^^ 

inspection  of  Figs.  1-3  shows  that  sub-unit  log  k^^  -  log  slopes  (i.e. 

the  fractional  power  law  k^^^  a  r'*,  where  0  <  x  <  1)  might  often  be 

anticipated  experimentally.  The  average  value  of  x  should  depend  both  on 

the  magnitude  of  and  the  range  of  values  accessed  (cf.  ref.  6b). 

The  presence  of  solvent  inertial  effects,  anticipated  for  low-friction 

solvents,  may  also  act  to  depress  these  slopes  somewhat,  especially  since 

Is  anticipated  to  vary  with  the  solvent  in  a  qualitatively  similar 

22 

fashion,  albeit  to  a  markedly  smaller  extent  than  r'^. 

Interestingly,  log  k^^^  -  log  t"^  (or  related)  plots  having  similarly 

fractional  slopes  have  been  predicted  on  the  basis  of  other  considerations, 

such  as  the  influence  of  reactant  vibrational  distortions®  and  to  the 

22 

presence  of  additional  high-frequency  solvent  relaxation.  The 

experimental  observation  of  such  a  "fractional  power- law"  rate-r*^ 
dependence  in  at  least  one  case^^^  has  been  attributed  to  the  former 

O 

factor.  It  would,  however,  seem  at  least  as  likely  that  this  observation 
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is  actually  due  to  the  friction-dependent  emergence  of  nonadiabatic 
effects,  as  considered  here.  In  any  case,  it  is  evident  that  such 
fractional  log  -  log  dependencies  should  often  be  expected.  This 
is  especially  true  given  the  likelihood  that  many  outer-sphere  electron- 
transfer  reactions  involve  relatively  weak  donor-acceptor  orbital  overlap, 
such  that  H°2  <  0.1  to  0.2  kcal  mol*^.^^’^^ 

In  order  to  investigate  more  closely  the  role  of  the 
functionality  upon  the  predicted  curves  it  is  also  instructive  to 

examine  the  dependence  of  k^^^  upon  r  for  varying  degrees  of  solvent 
friction  and  orbital  overlap.  Figures  4A  and  B  contain  plots  of  k^^ 

ratioed  to  the  value  when  r  -  r^^  (i.e.  reactant  spheres  in  contact), 
k^^/k°j^,  against  the  radial  separation  distance,  (r  -  r^)  .  The  rate 

constants  were  calculated  as  in  Fig.  2B,  thereby  accounting  both  for 
barrier-top  roundedness  and  solvent  inertial  effects.  The  four  traces 
shown  in  both  Figs.  4A  and  B  refer  to  values  of  0.02,  0.2,  0.5  and 
1.0  kcal  mol'^,  as  indicated.  The  values  in  Figs.  4A  and  B  are  held 
constant  at  5  x  10'®  s  and  5  x  10'^^  s,  respectively.  Although  somewhat 

arbitrary,  these  values  are  chosen  so  to  correspond  to  conditions  where 

the  reaction  remains  entirely  adiabatic  (at  r  -  r^^)  over  the  entire  range 
of  values  selected  (Fig.  4a) ,  and  where  the  reaction  switches  from 

essentially  nonadiabatic  to  adiabatic  behavior  as  Hjj  is  increased 

(Fig.  4B) .  Differences  in  these  ^  profiles  as  is  varied 

reflect  variations  in  the  relative  spatial- dependent  contributions  to  the 
overall  integration  yielding  k^^^  [Eq.  (1)). 

In  the  adiabatic  case  (Fig.  4A) ,  the  form  of  this  ‘ ^ 

dependence  is  seen  to  depend  only  to  a  small  extent  upon  The  slightly 
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more  pronounced  decrease  in  increasing  r  for  larger  values 

is  due  primarily  to  the  concomitant  increases  in  AG  caused  by  the  spatial 
variations  in  Slightly  greater  alterations  in  the 

profiles  with  increasing  Hjj  are  seen  in  Fig.  4B,  in  that  the  shallowness 
of  the  descent  with  increasing  r  is  maximized  at  an  intermediate  H®^  value, 
0.2  kcal  mol*^.  This  nonmonotonic  behavior  is  due  to  the  competing 
influences  of  the  Hjj-r  dependence  upon  and  AG*;  only  the  latter  factor 
influences  the  profiles  in  Fig.  4A.  In  any  case,  no  dramatic  changes  in 
these  profiles  with  varying  H°2  (or  are  seen  in  either  Fig.  4A  or  B. 

This  relative  invariance  of  the  shape  of  the  -r  profiles  is 

responsible  for  the  approximate  applicability  of  the  simple  encounter 
preequilibrium  treatment  in  the  present  context,  since  then  the  effective 
value  of  Kp  will  not  vary  greatly  with  either  H°2  or  This  perhaps 

fortunate  situation  is  brought  about  in  part  by  the  increases  in  AG  with  r 
predicted  from  Eq.  (5)  and  especially  Eq .  (6).  In  order  to  illustrate  this 
point,  a  set  of  (k,^A°t)'^  curves,  corresponding  to  the  conditions  of 
Fig.  4A  except  that  AG^  is  taken  to  be  independent  of  r,  are  shovm  in 
Fig.  5.  In  contrast  to  Fig.  4A,  the  (k,^A°t)-r  profiles  in  Fig.  5  extend 
to  substantially  larger  r  values  and  are  much  more  dependent  upon  H°2 . 
These  differences  graphically  demonstrate  the  importance  of  the  anticipated 
spatial  dependencies  of  AG*  to  the  integral  in  Eq.  (1). 

As  already  noted,  Beretan  and  Onuchic  have  recently  considered  in 
detail  some  factors  that  influence  the  form  of  the  bimolecular  rate 
expression  and  its  dependence  on  solvent  friction  in  the  intermediate 
regime  between  the  adiabatic  and  nonadiabatic  limits. Included  in 
ref.  15  is  a  discussion  of  some  assumptions  embodied  in  the  well-known 
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spatial  integration  formula  for  (Eq.  (1)].^^^  In  essence,  their 
analysis  utilizes  a  relation  similar  to  Eq .  (23)  so  that  no  explicit 
consideration  of  barrier-shape  effects  upon  k“j_  are  included;  as  noted 
above  this  expression  is  valid  only  for  «  k^T.  The  present  results 
also  differ  from  those  in  ref.  15  in  that  the  numerically  important  effects 
of  spatial  variations  in  AG  are  not  contained  in  the  latter. 

Electrochemical  Reactions 

Although  the  foregoing  analysis  applied  to  electrochemical  reactions 
proceeds  along  the  same  lines  as  for  bimolecular  processes,  the  form  of  the 
resulting  log  k^^^  -  log  plots  can  be  significantly  different.  This  is 
not  only  because  of  the  differing  form  of  the  spatial  k^^  integration 
[Eq.  (2)  versus  Eq.  (1)]  but  also  in  view  of  the  milder  dependence  of  AG* 
upon  the  donor-acceptor  separation  in  the  electrochemical  case  (vide 
supra) . 

Figure  6  shows  typical  plots  of  log  k'^^  versus  log  obtained  for 
conditions  corresponding  to  those  in  Fig.  3,  but  where  the  AG  -r  dependence 
is  given  by  Eq .  (6a)  and  the  k^^'^  integration  utilizes  Eq.  (2).  As 
before,  the  solid  and  dotted  curves  refer  to  the  presence  and  absence, 
respectively,  of  the  effects  of  barrier-top  roundedness  and  solvent 

,  with  AG^  for 

infinite  surface-reactant  separation  (i.e.  r^  -♦  <*>)  set  at  5  kcal  mol'^. 

11c 

(These  values  are  typical  for  some  experimental  systems.  )  The 
corresponding  dashed  curves  in  Fig.  6  were  obtained  from  the  encounter 
preequilibrium  treatment  as  in  Fig.  3,  but  with  taken  as  0.5  A  (cf  refs, 
lib,  13c).  As  before,  the  lower  and  upper  trio  of  curves  (displaced  as 


such  at  large  refer  to  H°2  values  (i.e.  Hjj  at  r^  -  8  A)  of  0.05  and 

0.5  kcal  mol'^. 
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Close  comparison  of  Fig.  6  with  Fig.  3  shows  that  the  r^-dependent 
log  -  log  t'^  slopes  for  the  corresponding  electrochemical  and 
homogeneous  curves  calculated  using  k^^^  spatial  integration  are  very 
similar  (within  ca  5%),  being  smaller  than  the  slopes  for  the  companion 
"encounter  preequilibrium"  curves  co  a  comparable  extent.  This  result 
might  at  first  sight  appear  surprising,  since  the  presumed  increases  in  AG* 
with  increasing  r  are  milder  in  the  electrochemical  case,  yielding  k^j^-r 
profiles  which  extend  to  larger  r  and  exhibit  a  greater  sensitivity  to 
than  for  the  homogeneous -phase  case  (vide  supra).  Offsetting  this  factor 
is  the  involvement  of  linear  spatial  coordinates  in  the  electrochemical 
k^j^-r  integration  rather  than  spherical  coordinates  as  in  the  homogeneous - 
phase  case.  In  comparison  with  the  latter,  the  former  function  is  weighted 
more  heavily  towards  rate  contributions  for  smaller  donor-acceptor 
separations . 

Overall,  then,  similar  log  k^^^  -  log  r‘^  dependencies  might  often  be 

anticipated  for  electrochemical  and  homogeneous -phase  exchange  processes. 

Differences  can  nevertheless  be  expected,  especially  since  the  values  of 

Hj2  and  their  functional  dependence  upon  r  may  be  quite  dissimilar  in  the 
35  36 

two  environments,  ’  Another  complication  for  electrochemical  process  is 

that  the  dependence  of  AG^  upon  r^  may  be  quite  different  to  that  predicted 

37 

by  the  simple  Imaging  expression  Eq.  (6a). 

Concluding  Remarks 

It  should  be  recognized  that  the  assumptions  embodied  in  the  foregoing 


analysis  impart  to  it  less  than  general  applicability.  In  particular,  the 

barrier-crossing  frequency  may  often  be  influenced  importantly  by  inner- 

shell  as  well  as  solvent  reorganization  dynamics.  The  analysis  will  also 

need  to  be  modified  for  reactions  having  large  driving  forces,  and  where 

electrostatic  work  terms  are  present  since  they  can  alter  substantially  the 

form  of  k^j^(r).  These  complications  will  nevertheless  be  absent  for 

homogeneous  self -exchange  processes  where  one  of  the  reacting  partners  is 

uncharged.  This  circumstance  prevails  in  the  experimental  companion  paper 

lie 

to  the  present  work.  °  Even  in  the  absence  of  such  work  terms,  however, 

the  assumptions  that  gjj^  is  independent  of  r,  and  that  Hj2  remains  constant 

when  the  solvent  is  varied  may  prove  to  be  significant  weaknesses  once 

specific  solvation  effects  upon  the  approach  of  the  reaction  partners  are 

considered.  The  desirability  of  including  the  effects  of  reactant 

1  6  38 

orientation  upon  the  distance  profiles  is  also  clear.  ’ 

On  the  other  hand,  the  predicted  rate-friction  dependencies  are  not 

particularly  sensitive  to  most  details  of  the  kinetic  model  for  k^^^,  such 

as  the  form  of  used,  the  r-dependence  of  Hj2.  and  the  modifications  to 

k*j^  caused  by  barrier-top  shape  effects. 

The  effect  of  non-Debye  relaxation,  which  appear  to  exert  marked 

influences  upon  k^^^  in  a  number  of  solvents ,  can  be  incorporated 

22 

into  the  numerical  treatment.  However,  given  the  magnitude  of  this 

2  2 

effect  together  with  complications  in  its  physical  description  the 
analysis  of  combined  electronic  coupling- solvent  friction  factors  is 
perhaps  most  prudently  restricted  to  Debye-like  media.  Deviations  from  the 
dielectric  continuum  approximation,  manifested  most  simply  in  spatial 
variations  in  r^,  can  also  be  considered,  although  on  the  basis  of  a  mean 
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spherical  approximation  (MSA)  treatment^  the  influence  upon  can  be 

39 

relatively  small. 

Given  that  most  outer-sphere  electron- transfer  processes  are  likely  to 
involve  relatively  weak  or  moderate  degrees  of  electronic  coupling, 
however,  the  present  treatment  should  provide  not  only  a  useful  means  of 
describing  the  effects  of  varying  the  solvent  friction,  but  also  can  yield 
valuable  Information  on  the  magnitude  of  itself.  Just  such  an 

application  from  our  laboratory,  to  the  self-exchange  kinetics  of  various 

Up 

cobalt  and  iron  metallocenes,  is  presented  elsewhere.  ^ 
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Appendix  A 

As  noted  in  the  text,  we  utilize  the  treatment  for  adiabatic  barrier 

24 

crossing  described  by  Okuyama  and  Oxtoby,  involving  the  concept  of  a 
position-dependent  diffusion  coefficient,  D^.  This  requires  the 

calculation  of  a  pair  of  diffusion  coefficients,  and  Dj^,  which  describe 
the  passage  of  the  system  sequentially  up  the  reactant  well  (for  E  < 
and  over  the  barrier  top  (E^  <  E  <  E^),  respectively.  Outlined  here  is  the 
derivation  of  Eqs .  (15)  and  (16),  which  express  and  in  terms  of  the 
effective  friction  and  t  respectively  [Eq.  (14)].  The  problem 

involves  the  calculation  of  and  for  two  Brownian  oscillators  with 

natural  frequencies  and  acted  on  by  external  forces  -dV^/dE  and 

-dVjj/dE,  respectively.  We  use  the  results  of  refs.  24  and  25  derived  on 
the  basis  of  the  Generalized  Langevin  equation  (GLE)  for  Brownian 
oscillators  and,  more  specifically,  the  Generalized  Smoluchowski  equation 
for  the  corresponding  GLE. 

The  diffusion  coefficient  Dp  is  defined  as  the  long  time  limit  of 
Dp(t),  expressed  in  terms  of  the  position  autocorrelation  function  ;^p(t)  - 
<E(t)  E(o)>p/<[E(o)]2>p: 

Dp  -  lim  Dp(t)  -  lim  (^mj.w2)'^[+  d  In  Xp(t)/dt]  -  (/3mgr?*“)'^  (A.l) 

t  -♦  «o  t  -♦  ® 

where  the  plus  and  minus  signs  are  approximate  for  p  -  b  and  p  -  w, 
respectively. 

Since  we  are  concerned  here  only  with  Debye  solvents,  the  friction  rj 
is  independent  of  frequency  and  equal  to  (e.g.  ref.  4),  whereby  x„(t) 

and  determined  by  the  inverse  Laplace  transforms: 

X„(t)  -  {(s  +  ri)/{s^  -I-  »js  +  w^))  (A. 2) 

>rb(t)  -  |(S  +  f7)/s*  +  TjS  ■  w^)) 


(A. 3) 
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The  solutions  Xp(t)  can  be  written  in  the  general  form: 

XpCt)  -  A^p  exp(Sipt)  +  Ajp  expCSjpt)  (A. 4) 

where  S^p  and  Sjp  are  the  roots  of  the  quadratic  equations  in  the 
denominators  of  Eqs.  (A. 2)  and  (A. 3)  for  p  -  w  and  p  -  b,  respectively,  and 

Alp  -  (Sjp  +  »j)/(S,p  -  S  ) 

(A. 5) 

A2p  -  (S^p  +  *j)/(S2p  -  S,p) 

The  correlation  function  Xp(t)  for  the  "instability  region"  Vjj(q)  is 
determined  by  the  positive  real  root^^’^^  -  [(t)^  +  4w^)^  -  (r])]/2. 

Inserting  Xb(t)  (  -  A^^,  exp(Sjj,t)]  into  Eq.  (A.l)  yields 
leading  to  Eq.  (16). 

For  the  overdamped  and  critically  damped  friction  regimes  (r;  >  2w^)  , 

24a 

X„(t)  is  determined  by  the  two  negative  roots  S2„  <  <  0  and  its 

long-time  behavior  is  determined  by  the  larger  root  i.e.  only  by  the 

first  term  in  Eq.  (A. 4).  Inserting  again  X„(t)  into  Eq.  (A.l)  yields 
^ Eq.  (15a). 

For  the  underdamped  regime  (rj  <  ,  we  obtain  two  conjugate  roots 

Siw  “  sj,  -  +  ib^:  also  A^^  -  aJ,  -  x  -e  iy.  Thus:^^^ 

X„(t)  -  2(x^  +  y^)*^  exp(a^t)  cos((p  +  b^t)  (A.  6) 

and 

I>„(t)  -  (^mjw2)'^[-a„  +  +  b^t)J  (A.  7) 

with  fp  -  arctg  (y/x) .  The  long  time  behavior  of  D^(t)  can  be  considered  to 

be  dominated  by  -  due  to  the  oscillatory  character  of  the  tg- 

24a 

term.  Making  this  approximation,  we  obtain  Eq.  (15b). 

It  is  also  of  interest  to  examine  the  forms  of  Eqs.  (15)  and  (16)  in 
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the  limiting  case  where  motion  in  both  the  well  and  the  barrier-top  regions 
are  overdamped,  i.e.  when  t)  »  and  17  »  respectively.  In  this  case, 
the  square  root  terms  in  both  Eq.  (15a)  and  (16)  can  be  approximated  as 
(1±  x^)^  ~  1  -  0.5x^,  where  x  -  2u^/t]  or  2u^/r} ,  respectively,  reduces 
these  relations  simply  to  -  r).  This  limiting  case  corresponds 
to  the  Kramers  diffusion  treatment  generalized  to  allow  for  barrier- top 
curvature,^  which  we  utilized  for  the  numerical  calculations  in  ref.  12. 


Appendix  B 

Considered  here  is  a  justification  of  the  form  of  [Eq.  (20)] 

employed  in  the  above  analysis. 

We  can  generally  express  in  the  form  [cf  Eq.  (21)] 


n 


et 


1  +  (Kl/KO 


-  iTir  ‘'n  > 


(B.l) 


From  Eqs.  (24)  and  (25),  the  "adiabaticity  criterion"  O  (  - 
«  kjT  [i.e.  where  Eq .  (B.l)  reduces  to  Eq .  (23)]: 

n  -  ^-1  (B.2) 


since  [cf  Eqs.  (24), (25)]: 

i/„  -  (B.3) 

Combining  Eq.  (B.l)  with  Eq.  (3)  yields 


K,^-n/(l-Hn)  (B.4) 

19 

The  well-known  form  of  [Eq.  (20a)],  suggests  the  use  of  the  related 

quantity  y^: 
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7,  -  (E^  V)""  (B.5) 

where  is  replaced  by  2it  in  order  to  allow  for  the  presence  of  solvent 
friction  (cf  ref.  40).  From  Eqs.  (B.2)  and  B.5)  we  obtain  n  -  4jr®^^7j; 
from  Eq.  (B.4)  we  can  express  as 

The  form  of  Eq.  (B.6)  suggests  intuitively  the  corresponding  expression 

[Eq.  (20)]  for  given  in  the  text. 

The  form  of  Eq.  (B.5),  and  its  TST  analog  Eq.  (20),  although 

approximate,  provides  a  useful  means  of  interpolating  between  the  adiabatic 

and  nonadiabatic  limits.  Table  B1  contains  a  numerical  comparison  between 

values  of  evaluated  for  typical  values  of  and  Hjj  (at  AG  -  6  kcal 

mol'^)  using  the  approximate  form  Eq.  (20)  and  the  exact  expression 

Eq.  (19).  Also  included  in  Table  B1  are  corresponding  values  obtained 

13a  14a 

using  the  oft-encountered  multiple-crossing  expression,  '  Eq .  (21). 
Inspection  of  Table  B.l  shows  that  Eq.  (20)  provides  a  generally  better 

approximation  than  Eq.  (21)  to  the  exact  expression  Eq .  (19),  although 

the  differences  are  not  large.  The  preference  for  such  a  Fade  form  for 
as  in  Eq.  (20),  has  also  been  noted  elsewhere. 

This  numerical  success  of  Eq.  (20)  also  underscores  the  general 
usefulness  of  Eq.  (B.6).  Thus  in  the  presence  of  friction,  Eq.  (B.6) 
becomes  equivalent  to  the  transmission  coefficient  extracted  by  inserting 
Eqs.  (23)  and  (B.3)  for  and  respectively,  into  Eq.  (3).  In  the 

absence  of  friction,  Eq.  (B.6)  reduces  to  Eq.  (20).  Equation  (B.6) 
therefore  provides  a  useful  approximate  means  of  accounting  in  a  combined 
manner  for  the  effects  of  friction  and  reaction  nonadiabaticity .  A  related 


approach,  involving  the  use  of  Eq.  (21)  with  wJI-k  replaced  by  was 

utilized  in  our  earlier  numerical  calculations  of  rate-friction 

^  ^  .  12.41 

dependencies . 


TABLE  Bl.  Numerical  Comparisons  of  the  Electronic  Transmission  Coefficient 
Jc  Evaluated  using  Different  Expressions  for  6.0  kcal  mol'^  barrier  at 


«o/2 

kcal 

n ,  s~ 

mol'^ 

1 

0.6 

10^® 

0.2 

10“ 

0.2 

10“ 

0.1 

10“ 

0.1 

19“ 

0.05 

10“ 

0.01 

10® 

^.1 

[Eq. 

(19)] 

0.504 

0.103 

0.540 

0.029 

0.219 

0.945 

0.785 

(Eq. 

(20)] 

0.520 

0.110 

0.553 

0.030 

0.236 

0.885 

0.756 

Si 

[Eq. 

(21)] 

0.590 

0.113 

0.631 

0.030 

0.250 

0.989 

0.881 
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Logarithmic  plots  of  the  calculated  rate  constant  for  electron  self 
exchange,  (M"^  s'^)  ,  against  the  inverse  longitudinal  solvent 

relaxation  time,  .  for  various  values  of  the  electronic  coupling  matrix 
element,  Values  of  calculated  using  the  simplified  expression  for 

(23)]  along  with  Eqs .  (1),  (4) -(6).  The  values  of  reactant 

contact,  H°2  1  are  0.02,  0.05,  0.1,  0.2,  0.5,  and  1.0  kcal  mol'^, 

corresponding  to  progressively  increasing  k^^^  values.  Other  parameters 
used  in  these  calculations  are  AG*  -  5.0  kcal  mol"^,  a  -  1.25  A‘^, 

r  -  2a  -  8  A,  and  T  -  298  K. 

Fig.  2 

As  in  Fig.  1,  except  that  k^^^  is  calculated  using  the  more  complete 
expression  [Eq.  (27)],  with  k^®^  given  by  Eq.  (18),  from  Eq .  (20),  and 
k“j^  from  Eqs.  (8) -(16).  Values  of  0)^/2n  taken  to  be  1  x  10^^  s'^  and  5  x 
10^^  s”^  in  Figs.  2A  and  2B,  respectively. 

£-iiL-3 

Comparisons  between  log  k^^^  -  log  plots  calculated  by  using 

spatial  integration  of  k^^^  (solid,  dotted  curves)  with  those  obtained  from 
the  simplified  encounter  preequilibrium  treatment  (dashed  curves) .  The 
lower  and  upper  sets  of  curves  refer  to  values  of  0.05  and  0.5  kcal 
mol*^,  respectively.  The  solid  curves  were  calculated  exactly  as  in 
Fig.  2B,  and  the  dotted  curves  as  in  Fig.  1.  The  dashed  curves  were 
calculated  using  k^^^  values  obtained  as  for  the  corresponding  solid  curves, 
but  by  inserting  k,^  for  r  -  (  -  8  A)  into  the  encounter  preequilibrium 

expression  [Eq.  (7)],  with  taken  as  0.25  (see  text). 
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Fig-  ^ 

Plots  of  normalized  to  the  value  with  reactants  in  contact, 

against  the  radial  separation  distance,  (r  -  r^^).  Values  of  k^^^ 
calculated  as  noted  in  caption  to  Fig.  2B.  Each  trace  refers  to  the 
different  values,  as  indicated  (kcal  mol'^).  Values  of  fixed  at 

5  X  10'*  s  and  5  x  10'^*  s  in  Figs.  4A  and  B,  respectively. 

Fig  ■  5 

★ 

As  in  Fig.  4A,  but  calculated  with  AG^  value  held  independent  of  r. 

Comparisons  between  log  k'^^  -  log  plots  for  electrochemical 

exchange  calculated  by  using  spatial  integration  of  k^^^  (solid,  dotted 
curves)  with  those  obtained  using  the  encounter  preequilibrium  treatment 
(dashed  curves)  (cf  Fig.  3  for  homogeneous  bimolecular  case).  The  lower 
and  upper  sets  refer  to  Hjj  values  of  0.05  and  0.5  kcal  mol'^, 
respectively.  The  solid  curves  take  into  account  the  effects  of  barrier 
top  roundedness  and  the  inertial  limit  upon  k^^^  using  Eq.  (27),  along  with 
Eqs,  (8-16),  (18),  and  (20).  The  parameters  used  are: 

AG*  -  5.0  kcal  mol'^  [for  r^  -  *,  Eq.  (6a)],  a  -  1.25  A'^,  r^^  -  2a  -  8  A, 
Wj,/2»r  —  5  X  10^^  s'^,  and  T  —  298  K.  The  dotted  curves  were  calculated 
using  Eq.  (23)  in  place  of  Eq.  (27).  The  dashed  curves  were  calculated 
using  k^j^  values  obtained  as  for  the  corresponding  solid  curves,  but  by 
inserting  k^^^  for  r,  “  ^0  (  •-  8  A)  into  the  encounter  preequilibrium 

expression  [Eq.  (7)],  with  Kp  taken  as  0.5  A  (see  text). 
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